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STAR Abstract 


A preliminary analysis is performed of deployment schemes and 
resulting impacted errors of orbiter-launched penetrators at bodies 
without atmospheres. Penetrators are missile shaped objects designed 
to successfully implant subsurface experiments on planetary bodies at 
impact speeds of about 150 m/sec. Results are presented for lunar and 
Mercury penetrator missions. In both cases it is shown that the 
deployment mass and impact errors are minimized if the penetrators 
are dropped at low periapse orbit altitudes. Closed-loop guidance is 
required to monitor errors in the large penetrator retro braking 
maneuvers in order to subsequently align the penetrator for near-zero 
angle of attack impact. 
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FOREWORD 


This study was performed between November 1974 and February 
1975 as part of the task schedule completed by Science Applications, Inc. 
for the Planetary Programs Division of OSS /NASA under Contract No. 
NASW-2613. The results are intended to assist NASA in advanced planning 
of planetary surface missions. Specifically, the objective of the report is 
to provide a preliminary evaluation of the feasibility of penetrator deploy- 
ment at targets without atmospheres. If identified control problems are 
not unreasonably complex, further detailed analysis can be supported to 
bring the gasless planet penetrator design to a level of design definition 
compatible with planning requirements. 

The author expresses his appreciation to E. Reese of Sandia 
Laboratories, to F. Jordan and C. Kohlhase of JPL for valuable technical 
inputs, and to A. Friedlander and J. Niehoff of SAI for many stimulating 
conversations and suggestions. 



SUMMARY 


Penetrators are missile shaped objects designed to implant electronic 

instrumentation in a wide variety of surface materials with a nominal 

impact speed around 150 m/sec. They have been used successfully in 

many terrestrial applications over the past decade. Recently they have 
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also been proposed for post -Viking /7 5 Mars exploration. ’ 7 The most 
significant advantage of planetary penetrators is that they avoid the high 
cost of soft landers without imposing the extreme impact conditions of 
hard surface landers on the payload. An initial favorable response by the 
science community to the exploration potential of Mars penetrators has 
prompted an interest in the application of this concept to in situ subsurface 
studies of other terrestrial bodies and planetary satellites. Unlike Mars 
many of these objects do not have atmospheres. A first order feasibility 
question has thus arisen: "Can penetrators be successfully guided to the 
required near-zero angle -of- attack impact conditions in the absence of an 
atmosphere ?" A preliminary answer to this question is the purpose of this 
report. 

The scope of the analysis includes two potential targets, i. e. the 
moon and Mercury, involves several different penetrator deployment modes, 
and focuses on impact errors arising from open-loop and closed-loop 
deployment control systems. Successful penetrator implacement requires: 
1) that the impact speed be controlled, nominally to 150 m/sec; 2) that the 
penetrator angle of attack, measured between the longitudinal axis and 
velocity vector, be in the range 0° - 11° at impact; and 3) that the impact 
flight path angle be within 15° of vertical. It is the errors in these 
terminal conditions that are the principal concern of this study. 
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The best mode of penetrator deployment identified uses an orbiting 
spacecraft as a penetrator launch platform. This mode, labeled the 
Intermediate Ellipse Transfer (IET) Mode, is depicted in Table SI. Prior 
to deployment the orbiter is first placed in an elliptical low-periapse 
altitude orbit. The penetrator is launched at periapse with a retro motor 
which kills its orbital velocity. It is then pitched over and allowed to free- 
fall to the surface. The deployment characteristics of the IET Mode are 
also summarized in Table SI. The initial orbit is circular at the moon and 
elliptical at Mercury. A free-fall impact velocity of 150 m/sec means that 
the periapse altitudes of these orbits must be lowered to 7 km and 3 km at 
the moon and Mercury, respectively, prior to penetrator deployment. 

This should not be a problem at the moon, but at Mercury a combination of 
several orbital maneuvers, onboard radar altimetry, and solar perturbation 
control will be necessary to achieve the very low altitude of 3 km. Also it 
will be undesirable to leave the orbiter in this orbit for more than several 
revolutions due to the impact hazard without continuous control. The 
penetrator retro AV requirements to kill the orbital periapse velocity are 
a nominal 1700 m/sec at the moon, but over 4 km/sec at Mercury. 

Assuming a penetrator impact mass of 31 kg (i. e. the Mars design), 10 kg 
for the attitude control system, and a single stage solid motor retro 
system, the total deployed mass of each penetrator is 81 kg at the moon 
and 251 kg at Mercury. These values can be compared with a deployed 
mass of only 50 kg at Mars where atmospheric braking is used to slow the 
penetrators. 

Key results of the deployment error analysis are summarized in 
Table S2. The errors in DSN tracking of the orbiter 's state at deployment 
are small and have little effect on any of the impact conditions except 
impact location. The primary error source for impact velocity and angle- 
of-attack errors is penetrator retro execution errors. The execution 
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Table SI 


REFERENCE DEPLOYMENT MODE 


• INTERMEDIATE ELLIPSE TRANSFER (IET) DEPLOYMENT SCHEMATIC 



• IET DEPLOYMENT CHARACTERISTICS 

Moon Mercury 

Initial orbit periapse altitude (km) 100 600 

Initial orbit eccentricity 0.0 0.8 

Penetrator deployment periapse altitude (km) 7.0 3.0 

Penetrator retro impulse (m/sec) 1698 4065 

Penetrator impact velocity (m/sec) 150 150 

Penetrator impact mass (kg) 31 31 

Penetrator deployment mass (kg) 81 252 


a. Just prior to single stage solid retro. 
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Table S2 


PENETRATOR DEPLOYMENT ERROR SUMMARY 3- 

Moon Mercury 


• DSN TRACKING ERRORS AT DEPLOYMENT 


Altitude (m) 

100 

150 

Velocity (m/sec) 

10 

24 

RETRO EXECUTION ERRORS (m/sec) 

Radial (x) 

30 

71 

In Path (y) 

26 

61 

Cross Path (z) 

30 

71 

OPEN LOOP IMPACT ERRORS 

Speed (m/sec) 

8 

41 

Angle of Attack (deg) 

15 

36 

Miss Distance (km) 

15 

21 

CLOSED-LOOP IMPACT ERRORS 

Speed (m/sec) 

8 

41 

Angle of Attack (deg) 

2 

2 

Miss Distance (km) 

15 

21 


a. 3 a errors of IET deployment mode 
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errors shown in the table are scaled to the magnitude of the impulse 
assuming 3a pointing errors of 1.0° and 3a magnitude errors of 1.5%. 

The affect of these errors on impact conditions are shown as open-loop 
impact errors. The critical errors are in impact angles of attack which 
are dominated by errors in the terminal flight path angle. With a maximum 
acceptable impact angle of attack of 11° required to successfully penetrate 
even very soft soils, it is readily seen that the open-loop control mode is 
unsatisfactory having 3a values of 15 degrees at the moon and 36 degrees 
at Mercury. Adding an accelerometer triad to the penetrator to monitor 
the retro burn errors easily reduces the angle of attack errors to very 
small values (2°) as can be seen by the tabulated closed-loop error 
summary. It should be noted, however, that nothing is done in the closed- 
loop mode to correct the execution errors, the attitude control system just 
accomodates them. Hence, the impact flight path angle, at Mercury in 
particular, may still be larger than the 15° off-vertical limit desired by 
some of the penetrator experiments, e.g. seismometers. 

As an overall conclusion to this analysis, the deployment of lunar 
penetrators appears to pose no unreasonable performance or control 
requirements. Conversely, the low deployment altitudes, the large retro 
mass, and large retro execution errors all raise feasibility questions for 
a Mercury penetrator mission. More detailed analysis will be required 
to resolve these issues and is recommended. 
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INTRODUCTION 


The use of subsurface penetrators in a variety of terrestrial 
scientific applications has been carried out over the past decade. These 
applications have included polar ice depth measurements and ocean floor 
sediment studies, along with more conventional impact studies into a 
variety of surface materials. Recent studies have considered the applica- 
tion of penetrators to the subsurface exploration of Mars as part of a post- 
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Viking/75 Mars scientific mission. ’ ’ 

The rationale for such missions is presented in references 1 and 2. 
Penetrators could return data on upper crustal heat flow, seismic 
measurements, physical properties measurements, geochemistry, and 
other in situ surface characteristics depending on the specific application. 
Successful penetrator deployment requires penetration to at least a one- 
meter depth in the hardest materials with maximum peak decelerations 
held to less than 2000 earth g T s for payload integrity. Penetrators were 
designed to meet these constraints for a wide variety of surface material. 
The Mars penetrator has a detachable afterbody which contains an 
antenna that remains on the surface and serves as a relay link between 
the buried forebody and the orbiting parent spacecraft. The forebody 
penetrates to a depth between 1 and 15 meters, depending on the surface 
materials. Table 1 gives a brief summary of the Mars penetrator as 
described in reference 1. To achieve the desired penetration depths 
requires an impact speed of about 150 m/s. The impact flight path angle 
should be greater than 75° to insure proper seismic measurements with 
fixed body detectors. Finally, the impact angle of attack (the angle 
between the longitudinal axis of the penetrator and the velocity vector at 
impact) is constrained to a small value. As the penetrator very quickly 
aligns itself with the velocity vector after impact, misalignment will give 
rise to large shear forces which could lead to penetrator failure. The 
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Table 1 


MARS PENETRATOR CHARACTERISTICS 


MASS SUMMARY 


at orbiter deployment 

50 

kg 

at impact 

31 

kg 

science/electronics budget 

7 

kg 

LENGTH 

140 

cm 

PRINCIPAL DIAMETER 

9 

cm 


IMPACT CONSTRAINTS 

speed 
direction 

angle of attack (rock) 

(loess) 

PEAK DECELERATIONS 

forebody (science) <1800 g for 5 ms 

afterbody (transmitter) <18000 g for 3 ms 


140-170 m/s 
vertical ±15° 
<3° 

< 11 ° 


DEPTH OF PENETRATION 


1-15 m 



magnitude of the allowable misalignment depends significantly on the 
impacted material. For rather unconsolidated surfaces such as loess, 
the alignment error could be as much as 11°, whereas for rocky surfaces 
3° would be a more reasonable upper bound. These penetrator impact 
constraints are also summarized in Table 1. 

For the Mars studies, the penetrator may be deployed from either 
a Pioneer or Mariner class spacecraft in orbit about Mars. The pene- 
trator and its aeroshell are inserted on an entry trajectory into the 
Martian atmosphere using a small solid rocket motor. Attitude control 
and orientation is provided by the parent spacecraft for this maneuver. 
Atmospheric deceleration and entry heating are controlled by the aero- 
shell design and the atmospheric entry angle. After the descent velocity 
becomes subsonic a secondary deceleration system, e. g. a parachute, 
is deployed for final terminal velocity and attitude control. 

The penetrator science return is probably higher for a body such 
as Mercury than it is at Mars, because no prior in situ subsurface 
measurements are likely to have been made. Initial seismic and heat 
flow measurement would greatly increase our knowledge of Mercury’s 
interior and its massive heavy core. For lunar applications, although 
additional measurements might not be as significant as the initial Apollo 
and Surveyor data, they would allow sampling of different lunar geologic 
features and permit a comparison of heat flow in the lunar uplands with 
the Apollo mare measurements. Seismic measurements could supple- 
ment or extend those currently being returned from the ALSEP stations. 
Hence, there exists a general interest in the broader application of the 
penetrator concept to solar system exploration. 
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The objective of this study is to examine the errors associated with 
penetrator deployment at bodies without atmospheres, with Mercury and 
the moon being emphasized in the analysis. For such bodies, the impact 
speed, flight path angle and attitude must be controlled through propulsive 
maneuvers and systems associated with the penetrator. Initially the 
penetrator deployment system is assumed to be "open loop”, i. e. no 
onboard measurements are taken and errors in impact conditions are 
calculated by propagating the errors resulting from, a) the orbit deter- 
mination process, b) maneuver execution, and c) attitude control system 
pointing, down to impact. It will be shown that an "open-loop” deploy- 
ment scheme cannot meet the impact constraint requirements presented 
in Table 1. Consequently, a "closed-loop" approach is also investigated 
using several different types of onboard measurements. This study is 
not intended to be a mission analysis report, but rather to ascertain the 
magnitude of the problem for penetrator deployment onto gasless bodies. 
Hence, the impact of orbital operations necessary or desirable for 
penetrator deployment on other mission objectives has not yet been 
examined. Subsequent sections of this report describe the penetrator 
deployment modes, error sources and error analysis, and the types of 
onboard measurements necessary to reduce terminal errors. 


/ 
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2 . 


PENETRATOR DEPLOYMENT MODES 


The penetrator impact constraints (see Table 1) require a near 
vertical approach to the surface with a planned terminal velocity of 
150 m/sec. These conditions imply a rectilinear free-fall descent from 
some altitude at which the penetrator is brought to rest. This rest 
altitude will depend on both the impact body mass and the desired impact 
speed. Rest altitudes and free-fall times are given in Table 2 for two 
impact speeds, 150 m/sec and 300 m/sec. The lower value of 
150 m/sec is the Mars penetrator design point and is considered a good 
compromise between impact deceleration and penetration for a wide 
range of impact materials. The higher value of 300 m/sec is probably 
an upper bound on impact speed. For most rocky surfaces (with soil 
constants 1 less than 1. 5) peak decelerations would exceed 2000 g's in 
the forebody at this speed, and for loose soils such as loess (with soil 
constants greater than 10) the depth of penetration would exceed 20 m. 
These values assume the same goemetric design as the Mars penetrator. 

The finite intrinsic rotational speed of the impact body surface 
will add a normal component of velocity to the free-fall impact speed 
with the result being only a nearly rectilinear rather than truly 
rectilinear terminal descent. For the two bodies of primary interest 
in this analysis, the moon and Mercury, their maximum surface rotational 
speeds are 4. 6 and 3. 0 m/sec, respectively, and occur, of course, along 
the equator. Using a minimum vertical impact speed of 150 m/sec, the 
maximum off-vertical impact angles are 1. 75° and 1. 15°, respectively. 

As these angles are small and would be even smaller for non -equatorial 
impact sites and higher impact velocities, the analysis and terminology 
of a true rectilinear descent is used in this report . 
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PENETRATOR FREE FALL CHARACTERISTICS AT ATMOSPHERELESS BODIES 
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Callisto 106.4 2500 1.14 


It is assumed that a parent spacecraft carries the penetrators into 
an initial orbit about the target body from which it can selectively deploy 
the penetrators and frequently communicate with them after they are 
implanted. Once in orbit a penetrator can be placed on a rectilinear 
descent ellipse (for a normal surface impact) at any time by killing its 
orbital velocity. The two extreme points of deployment considered in the 
analysis are at the orbit periapsis and apoapsis. The two modes 
associated with these orbit extremes are illustrated in Figure 1. The 
Intermediate Ellipse Transfer (IET) Mode first puts the spacecraft on a 
Hohmann transfer ellipse with a periapsis altitude equal to the desired 
rest altitude (see Table 2). The penetrator is deployed at periapsis with 
a AV which kills its orbital velocity and free falls to the surface. The 
second of the two modes, the Rectilinear Ellipse Transfer (RET) Mode, 
stops the penetrator at apoapsis from which it starts its free fall. As it 
nears the surface a second AV again stops the penetrator from which 
point (i. e. the rest altitude) it accelerates to impact at the desired speed. 

Deployment of penetrators from orbit with minimum impulse 
expenditure requires a low altitude near- circular orbit for either mode 
presented in Figure 1. Such an orbit is easily attained at the Moon, but 
at Mercury elliptical orbits are preferred due to the limited capacity of 
current propulsion systems to deliver useful payload into Mercury orbit. 
Reference initial spacecraft orbits for this analysis were selected with 
these considerations in mind; they are summarized at the beginning of 
Table 3. For lunar missions the initial orbit is circular with a low 
altitude of 100 km. For Mercury missions a 0. 8 eccentricity orbit with 
an initial periapse altitude of 600 km was chosen from reference 4. 


\ 
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FIG. 1. PENETRATOR DEPLOYMENT MODES 



The deployment impulse requirements at both the moon and Mercury 
are also summarized in Table 3 for both deployment modes presented in 
Figure 1, using both 150 and 300 m/sec impact speeds. From the data in 
the table it is readily seen that the IET Mode results in lower impulse 
requirements; the difference in impulse requirements between the two 
impact speeds is negligible. Also, the requirements at Mercury, due to 
the higher energy spacecraft orbit, are more than twice as large as the 
lunar cases. Variations from the nominal IET and RET Mode deployment 
schemes were briefly investigated with no significant difference found in 
the total impulse requirements. Included in the variations were 
1) performing AVA in the RET Mode 10° off apoapse, 2) employing a 10° 
off-normal terminal descent, and 3) reducing the initial Mercury apoapse 
altitude to 9700 km. The largest reduction in the deployment impulse, 
due to reducing the Mercury apoapse altitude, was less than 4% of the 
total requirement. 

Although the IET Mode readily produces the lowest impulse 
requirement it also creates the greatest risk to mission success since the 
spacecraft must come very close to the surface before releasing the 
penetrator, only 3 km at Mercury if an impact speed of 150 m/sec is 
required. In fact, drawing on results presented below in Section 4, the 
3aperiapse altitude dispersions resulting from the AVA maneuver are 
about 8.4 km. Consequently, there exists an unacceptable chance that 
the spacecraft would crash if it attempted to reduce the periapse 
altitude as required in one apoapse maneuver. This can be remedied by 
reducing the periapse altitude in steps, separated by periods of orbit 
tracking and radar altimeter measurements in the vicinity of periapse. 
Alternatively, at Mercury, solar perturbations might be used to advantage 
producing the same desired reduction in periapse altitude. 
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Table 3 


IMPULSIVE PENETRATOR DEPLOYMENT REQUIREMENTS 


• INITIAL SPACECRAFT ORBITS 


Moon Mercury 


Periapse Altitude (km) 100 

Apoapse Altitude (km) 100 

Eccentricity circular 

Period (hours) 2.0 


• DEPLOYMENT IMPULSE REQUIREMENTS (m/sec) 
o IET Mode - 150 m/sec impact 


AVA a (see Fig. 1) 21 

AVB 1698 

Total Impulse 1719 

o IET Mode - 300 m/sec Impact 

AVA a 16 

AVB 1683 

Total Impulse 1699 

o RET Mode - 150 m/sec Impact 

AVA 1632 

AVB 536 

Total Impulse 2168 


o RET Mode - 300 m/sec Impact 


600 

24880 

0.8 

22 


38 

4065 

4103 


37 

4057 

4094 


399 

4049 

4448 


AVA 

AVB 

Total Impulse 


1632 399 
471 4041 
2103 4440 


a. AVA maneuver is made by the spacecraft propulsion system for the 
IET Mode. 
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The magnitude and direction of periapse altitude changes resulting 
from solar perturbations depend upon the orbit geometry and orientation. 
To successfully reduce the periapse altitude to a very small value without 
impacting the surface requires that the orbit to orbit change become 
vanishingly small as the final value of periapse is approached. A 
numerical analysis of various Mercury orbit conditions was performed 
for mission launch opportunities in 1980, 1985, and 1988 to see if orbits 
with this characteristic exist. Results from the analysis are presented 
in Figures 2 and 3. From these figures it can be seen that the aim angle 
(measured clockwise in the impact plane from the T axis to the aim 
vector, B ) controls the direction of periapse change. Using these 
results, the periapse histories for two specific orbits are tabulated in 
Table 4 which have the desirable attributes of near zero change near the 
planet's surface. The data illustrate that orbits can be found which 
permit the consideration of very close approaches to Mercury's surface. 

In all the above discussion the principal concern has been the very 
low altitude requirement at Mercury required by the IET Mode in order 
to satisfy a 150 m/sec impact constraint. It should be noted from 
Table 2 that doubling the impact speed to 300 m/sec would raise the 
periapse altitude of the intermediate transfer ellipse from 3 km to a 
more comfortable 12 km, but would restrict the range of soils the 
penetrator could successfully impact. The evaluation of such a trade- 
off was not considered as part of this study. Suffice it to say that using 
solar perturbations as part of the periapse reduction strategy might 
reduce the impact hazard, but does add a new constraint on orbit 
orientation (particularly inclination and location of periapse) which will 
complicate the orbit selection process for the mission designer. 
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FIG. 2. AIM ANGLE EFFECT ON MERCURY ORBIT SOLAR 


PERTURBATIONS FOR THREE MISSION OPPORTUNITIES 
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FIG. 3. AIM ANGLE EFFECT ON MERCURY ORBIT SOLAR 

PERTURBATIONS FOR THREE ORBIT ECCENTRICITIES 








Table 4 


SOLAR PERTURBED PERIAPSE HISTORIES OF MERCURY ORBITS 


Launch Opportunities 
1980 1988 


• ORBIT PARAMETERS 


Period (hrs) 

22 

22 

Eccentricity 

0.8 

0.8 

Initial Periapse Altitude (km; 

138 

59 

Initial Epoch 

14 Apr 82 

17 Sep 90 

• PERIAPSE ALTITUDE HISTORIES (km) 

Orbit No. 0 

138.0 

59.0 

5 

132.6 

58.7 

10 

115.0 

55.7 

15 

85.4 

47.8 

20 

53.5 

36.3 

25 

28.2 

24.5 

30 

12.9 

15.1 

35 

5.9 

8.9 

40 

3.8 

5.6 

42 

3.6 

4.8 

44 

3.5 

4.3 

46 

3.5 

4.0 

48 

3.5 

3.8 

50 

3.3 

3.7 

52 

3.0 

3.7 

54 

2.4 

3.7 


a. Initial periapse altitude is reduced from the baseline altitude of 
600 km by a spacecraft apoapse maneuver. 



A preliminary estimate of the weight of the entire penetrator 
deployment system can be obtained by using the impact mass of the Mars 
penetrator design, i. e. 31 kg as presented in Table 1. To this should be 
added a 10 kg allowance for attitude control, sensors and control logic for 
guiding the penetrator between deployment and impact. For the IET Mode 
a minimum of only one retro motor is needed to kill the intermediate 
ellipse periapsis velocity. For the Mercury mission, however, a two- 
staged retro system will more effectively kill the large periapsis velocity. 
For the RET Mode at least two retro motors are needed if a simple solid 
propellant system is used. The first motor kills the orbital speed at 
apoapse and the second motor brakes the fall of the penetrator just prior 
to impact. A third motor might be used in the case of Mercury to split 
the large braking maneuver required. A summary of penetrator pre- 
deployment masses at the moon and Mercury is given in Table 5 for both 
deployment modes using the various staging strategies just discussed. A 
solid propellant specific impulse of 290 sec and a useful stage mass fraction 
of 90% has been assumed in generating the results which are for 
150 m/sec impacts. Compared with the 50 kg predeployment mass of the 
Mars penetrator, these data reflect much higher requirements - from 
almost two to nearly six times as much mass. Additional staging at 
Mercury does reduce the mass requirements by about 14%, but also 
complicates the maneuver sequence. On the basis of mass considerations 
alone, the IET Mode for lunar penetrators still appears tractable. At 
Mercury, however, considering the difficulty with placing mass into orbit 
with current propulsion systems, the large retro mass requirements may 
pose a significant limitation on the penetrator mission concept. 
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PENETRATOR PREDEPLOYMENT MASS (KG) SUMMARY 
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Based on Mars penetrator design; impact speed = 150 m/sec 

Assumes solid retro system; Isp = 290 sec; stage mass fraction = 0.9 and for multistage 
systems, 2% of the stage weight above the retro is allowed for interstage connection. 


3. 


ERROR SOURCES FOR PENETRATOR DEPLOYMENT 


As successful penetrator implantation must be achieved through 
impulsive maneuvers and active penetrator attitude control, it is 
imperative to understand how well terminal constraints can be met using 
active systems. Errors in the final conditions arise primarily from two 
sources: 1) uncertainty in our measurement of the basic variables, i. e. 
spacecraft state vector, planetary topography, etc. and 2) errors 
introduced by spacecraft systems, i. e. maneuver execution errors and 
attitude control errors. 

The orbit determination (O. D. ) error model used in this study was 
obtained from JPL, and represents the O. D. errors in DSN tracking of 
planetary orbiters after several initial orbits have been tracked and 
utilized to improve a priori estimates of gravitational parameters. 

Figure 4 defines the geometry used in this formulation. The x-y plane is 
in the plane of the orbit and z is parallel to the orbital angular momentum 
vector. The 3 ct uncertainties in the state vector are taken to be: 

Ay = . 006r, 

Az = . 006r, 

Ax ~ 10 -2 Ay, 

Ax = . 006 v cos y, 

Ay = . 006 v sin y, 

Az = . 006 v sin y. 

In these equations, r is the magnitude of the radius vector to the space- 
craft, v is the orbital speed of the spacecraft and y is its flight path angle. 
The errors are assumed to be uncorrelated. Note that at apoapse and 
periapse, where the flight path angle is zero, the uncertainty in the 
tangential component of velocity is very small, with essentially all of the 
uncertainty being in the radial component. This arises not from 
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uncertainty in the timing at which the spacecraft is at apsis, but rather 
from the uncertainty in the radius of the apsis location itself. The same 
formulation for O. D. errors is also used for circular orbits, i. e. for the 
lunar missions, and should be conservative for this case. Table 6 
summarizes the resulting O. D, uncertainties for the nominal orbits used 
in this study. 

The maneuver execution errors were formulated in terms of the 
variation in the magnitude and direction of the AV vector from nominal 
conditions, representing shut-off timing errors and the errors in the 
attitude control system in orientating and holding the spacecraft during 
the burns. The nominal 3a impulse error was taken to be 1. 5% of AV 
magnitude and orientation errors of 1.0° were adopted for both inplane 
and out-of-plane pointing. Table 7 gives the resulting velocity errors with 
both descent modes beginning from the nominal Mercury and lunar orbits, 
and assuming an impact speed of 150 m/sec. 

In a real application, the pointing errors will be related to the type 
of attitude control system used. A spin stabilized system is not well 
suited for penetrators since they are not dynamically stable when rotating 
about their roll axis. Furthermore, it would be difficult to pitch a spin- 
ning penetrator through 90°, particularly after the periapse impulse 
(AVB) of the IET Mode. Hence, three axis stabilization appears pre- 
ferable for penetrators based on initial consideration of the deployment 
requirements at atmosphereless bodies. 
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Table 6 

ORBIT DETERMINATION UNCERTAINTIES 8- 



Lunar 

Circular 

Orbit 

(100 km alt) 

Mercury e 

= 0. 8 Orbit 

Periapse 
(600 km alt) 

Apoapse 
(24880 km alt) 

Ax (km) 

0.1 

0.2 

1.6 

A y (km) 

11.0 

14.5 

164 

Az (km) 

11.0 

14. 5 

164 

Ax (m/sec) 

9.8 

24.4 

2.4 

Ay (m/sec) 

—0 

— 0 

~0 

Az (m/sec) 

—0 

— 0 

—0 


a. 3a values 


Table 7 

MANEUVER EXECUTION ERROR SUMMARY a 


Lunar Cases Mercury Cases 

IET Mode RET Mode IET Mode RET Mode 



AVA 

AVB 

AVA 

AVB 

AVA 

AVB 

AVA 

AVB 

Ax (m/sec) 

.4 

29.6 

28.5 

8.0 

.7 

70,9 

7.0 

61.0 

Ay (m/sec) 

.3 

25.5 

24.5 

9.3 

.6 

61.0 

6.0 

71.0 

Az (m/sec) 

.4 

29.6 

28. 5 

9.3 

.7 

70.9 

7.0 

71.0 


a. 150 m/sec impact speed 
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4. 


OPEN-LOOP ERROR ANALYSIS 


The parameters of interest (dependent variables) for penetrator 
deployment are speed, flight path angle at impact, and miss distance, 
defined as the surface range between the nominal impact point and the 
perturbed impact point. Also, errors in the altitude of the AVB point 
(see Fig. 1, pg. 8 ) for both the IET and RET Modes are of concern. For 
each variable of interest, q, a sensitivity matrix, giving the partial 
derivative of q with respect to each component of initial position and 
velocity is computed, i.e. 



M 

ai o 


aq 


i 




(1) 


The error sources described in Section 3 are utilized to construct 
covariance matrices of errors, a. These errors are assumed un- 
correlated, hence the covariance matrices, whose elements are the 3a 
error estimates, are diagonal. The variance in the dependent variable is 
then computed as 

a 2 (q)=S(q)AS(q) T . (2) 

Explicit expressions for all partial derivatives used are given in 
Appendix A. 

4. 1. Error Sensitivities 

Using the partial derivative expressions derived in 
Appendix A, numerical partial derivatives were computed for both the 
IET and RET Modes at the moon and Mercury. The results are presented 
in Table 8. The dependent variables are listed at the head of each column 
with the independent variables given in the left-hand columns. The partial 
derivatives of periapsis radius (RP) pertain to the IET Mode. Partials 
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Table 8 


PENETRATOR FLIGHT PARAMETER SENSITIVITIES TO 
DEPLOYMENT POSITION /VELOCITY 




Penetrator Flight Parameters 


Periapse 

Altitude, 

RP(km) a 

Braking 

Altitude, 

R(km)b 

Impact 

Velocity, 

v(m/s) 

Impact 
Angle, 
Y (deg) 

Miss 

Distance 

(km) 

• LUNAR CASES 






x (km) 

2.8 

1.1 

11 

-10- 6 

i 

o 

1 

o 

y n (km) 

— 10" 7 

0 

- 10- 7 

-10- 5 

1.0 

O 

z (km) 

— 10" 7 

0 

~ 10- 7 

— 10" 5 

1.0 

o 

x (m/sec) 

-10" 6 

0.4 

6. 7xl0 3 

-l<f 6 

— 10" 4 

o 

y (m/sec) 

4.2 

~ 10- 7 

6. 7xl0" 3 

0.4 

0.1 

U 

z o (m/sec) 

~ 10- 6 

c- 

i 

o 

T-t 

\ 

6. 7xl0 -3 

0.4 

0.1 

• MERCURY CASES 





x (km) 

0.2 

7.6 

25 

~ 10- 6 

o 

o 

t 

0 

y„(km) 

— 10 " 8 

0 

0.01 

~i<r 6 

1.0 

U 

z (km) 

-10" 8 

0 

0.01 

-10- 6 

1.0 

U 

x (m/sec) 

— 10" 6 

137.3 

6.7xl0“ 3 

-xo - 6 

-xo - 8 

U 

y (m/sec) 

14.6 

- 10- 6 

6. 7xl0~ 3 

0.4 

0.04 

z Q (m/sec) 

-10" 8 

~ 10- 6 

6. 7xl0 3 

0.4 

0. 04 


a. For the IET deployment mode 

b. For the RET deployment mode 
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of radius magnitude (R) pertain to the RET Mode. Both of these sets of 
partials are taken at the rest altitudes which result in a nominal impact 
velocity of 150 m/sec. The remaining sets of partials, for impact 
variables, pertain to both deployment modes. 

A number of important points are observed from the data in 
Table 8. These can be summarized as follows: 

a) The periapsis error for the IET Mode is the result of two initial 
errors at the AVA maneuver, 1) the error in the radius 
magnitude, x q , and 2) the error in tangetial speed, y Q ; 

b) With the RET Mode the error in the radius (and altitude) of the 
AVB maneuver is dominated by initial errors in, 1) the apoapsis 
radius, x q , and 2) the radial speed, x q , after the AVA 
maneuver. 

c) Impact speed errors result primarily from 1) errors in all 
components of velocity after the AVB maneuver, and 2) the error 
in the radius, x q , at the AVB point; 

d) The error in the impact flight path angle, y, of the rectilinear 
descent results primarily from errors in the two velocity 
components of the AVB maneuver, i. e. y Q and z q . Note that the 
error in impact flight path angle is equivalent to the impact 
angle of attack for an open loop controlled penetrator on a 
nominally rectilinear final descent; 

e) Errors in impact location, i. e. miss distance, result from the 

uncertainties in both the position and velocity components normal 

to the local radius at the completion of the AVB maneuver which 

include x and x , and y and y . 
o o’ J o J o 
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Taking all of these effects together, it becomes apparent that a 
final rectilinear descent, which is used in both modes, translates errors 
in every component of the state at completion of the AVB burn into some 
significant form of error in the penetrator impact conditions. 

4.2. Impact Error Magnitudes 

Insight into the magnitude of the final impact errors may be 
gleaned by combining the dominant partial derivatives of Table 8 with the 
magnitudes of the initial O.D. and execution errors cited in Tables 6 and 
7, respectively. An illustration of these error magnitudes is presented 
here for the IET Mode applied to the Mercury penetrator case. Just the 
largest terms are analyzed to simplify the approximation of final errors. 
Trends in the lunar case are completely analogous. 

a) Periapsis Magnitude 

The net effect of initial position errors is small. The 

component x q , which has the largest partial derivative, has the smallest 

error component while the largest errors, in y Q and z q , have negligible 

partial derivatives. By far the dominant term is due to the 0. 57 m/s error 

in y , arising from execution of AVA, coupled with J*L j*j_ = 14.6 sec. 

° dy n 

Hence 0 

Arp pi Ay = 8. 3 km. 

sy 0 

J o 

It should be noted that the above is the uncertainty in the 
magnitude of periapse radius, whereas for penetrator deployment it is the 
altitude which is critical. The altitude errors depend on both periapse 
errors and uncertainties in topography. However, the periapse error is 
so large that surface relative measurements are required as discussed 
in Section 5. 
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b) Impact Speed 

The impact speed error due to position errors at AVB is 
most sensitive to the altitude error, Ax, as seen from Table 8. Assuming, 
however, that the periapse altitude of 3 km has been approached slowly 
with continual DSN tracking, the error in Ax is only 

Ax ~ . 00006 R = 150 m, 

P 

which leads to an impact speed error of only about 4 m/sec. The execution 
errors at AVB are large, 61-71 m/sec as seen in Table 7, but at first 
glance the sensitivities given in Table 8 appear to be small enough to 
negate any significant error contribution. Unfortunately, the linear error 
analysis using partial derivatives with respect to initial velocity compo- 
nents are valid only for small initial errors. Using the sensitivities 
evaluated along the nominal trajectory with the large execution errors 
leads to a serious underestimate of the impact speed errors. 

For example, the sensitivity of impact speed, v, with 
respect to initial speed, v q , can be calculated by differentiating the energy 
equation as 

3V v q 

0 V o"v 

Evaluating this along the nominal trajectory, i. e. , v =1 m/s and 

dv -3 ® 

v = 150 m/s, gives 2 — = 7x10" . Using the components in Table 7, the 

„ ° v o 

3ct execution error, 6v , at AVB is 

’ o’ 

6v = 70. 9 2 + 61. 0 2 + 70. 9 2 ly/2 = 117 m/sec. 

o J 

The linear error theory gives for the impact speed error, 6v = 0. 8 m/s. 
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However, by applying the energy equation directly, the 
correct impact velocity is 190 m/s for a 117 m/sec error in initial speed, 
which amounts to a 40 m/sec impact speed error. To illustrate the non- 
linearity of this sensitivity, la execution errors lead to an impact speed 
error of only 5 m/sec. 

c) Impact Flight Path Angle 

This parameter is sensitive essentially only to the initial 
errors in y Q and z q . The dominant error here is the 61 m/sec error in 
y Q due to AVB, resulting in 

Ay« A-2X) Ay =23.4°. 

Wo) ° 

d) Miss Distance 

The error in the y and z components of position, along with 
the y error introduced by AVB, are the most significant, but the position 
errors, having larger partials, dominate the miss distance, D m . 

Hence 



where D is distance on the impact surface. 

The effects of all errors on the various parameters of 
penetrator deployment and impact are given in Table 9 for both the IET 
and RET Modes for the nominal Mercury case depicted in Figure 1, 
page 8 . The data were calculated from equation (2) using the computer 
program described in Appendix B except for the impact speed errors which 
were calculated as described above. The results presented include the 
total error in each parameter along with the contribution to the total error 
arising from each error source considered separately. 
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Table 9 

SUMMARY OF MERCURY PENETRATOR 3c ERRORS DUE TO O. D./MANEUVER UNCERTAINTIES 
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Path angle error equivalent to impact angle of attack 



Considering the IET Mode first, the error in the intermediate 
periapse radius (and altitude) of more than 8 km is due primarily to 
execution errors of the periapse lowering maneuver, AVA. Since a 
deployment periapse altitude of only 3 km is desired, this error is 
obviously much too large for an acceptably safe maneuver strategy. A 
combination of impulsive maneuvers, perhaps including solar perturbations, 
to reduce the periapse altitude, with an onboard spacecraft altimeter 
monitoring the periapse reduction, are probably required for this deploy- 
ment mode. Assuming this to be the case, the remaining three errors at 
impact, i.e. Av, Ay, and miss distance, result from O. D. errors at the 
AVB point and AVB execution errors. The only error of these three which 
is obviously unacceptable is the error in flight path angle, Ay, of 36°. 

Since, for open-loop control, the penetrator would be aligned after AVB to 
an attitude determined before AVB assuming a perfect maneuver, Ay would 
be equivalent to the actual angle of attack at impact. The maximum 
acceptable angle of attack in very soft material is about 11° (see Table 1, 
page 2 ). Hence, a successful impact would occur much less than even 
the la probability since there is an additional uncertainty in the hardness 
of the impact soil. 

The second half of Table 9 contains summary error data for the 
RET Mode. It is immediately apparent that all the total errors in this 
mode are larger than for the IET Mode. This is primarily due to the fact 
that the penetrator is deployed much further from the planet (at apoapsis) 
allowing errors to propagate from more sources and for a longer period of 
time. Both the errors in impact speed and flight path angle are unaccep- 
tably high. As a point of observation, the RET Mode is generally worse 
than the IET Mode from an open-loop control point of view. Finally, it 
should be noted that the miss distance from the targeted impact point is 
relatively small for both deployment modes. This small error is an 
inherent characteristic of rectilinear trajectories. 
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4.3. Detailed Open-Loop Error Analysis 


In an attempt to reduce the magnitude of penetrator 
deployment errors and to ascertain the sensitivity of these results to the 
various deployment input parameters, a number of off-nominal deployment 
cases were run for both the IET and RET Modes. Table 10a describes the 
conditions of each case and Table 10b and 10c give the resulting errors for 
lunar and Mercury missions, respectively. Two entries are given for each 
item within each case which pertain to the IET and RET Modes, respectively. 
For the IET Mode the deployment strategy includes the spacecraft lowering 
periapse to the desired altitude through a combination of impulsive 
maneuvers and solar perturbations and then releasing the penetrator. The 
errors are then computed based on tracking errors at the AVB point along 
with AVB execution errors. The comparable approach is not available for 
the RET Mode, since the spacecraft cannot deliver the penetrator to the 
AVB point, nor is it feasible to consider tracking, orbit determination, 
and a trajectory trim maneuver system for the penetrator after it has been 
released. Hence, errors in this mode are determined by propagating 
O. D. and execution errors at AVA to the AVB point and then adding AVB 
execution errors. The O. D. and execution error matrices at AVA have 
been propagated using the state transition matrix (see Appendix A). 

The Mercury cases (Table 10c) are discussed first. Moving AVA 
off apoapsis makes the final errors in periapse conditions significantly 
worse. Changing the size of the initial orbit results in a slight improve- 
ment in the IET Mode, and significantly reduces the terminal errors in the 
RET Mode, although not to an acceptably small level. Reducing the mag- 
nitude of the error sources by a factor of 3 shows final errors also 
reduced by a similar factor. Error levels of this reduced magnitude, i. e. 

0. 5% (3a) in A V magnitude and 0. 33° pointing accuracy, would be difficult 
to achieve in practice. Nevertheless the error levels, except for RP and 
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Table 10a 


DETAILED ERROR ANALYSIS SUMMARY - CASE DESCRIPTIONS 


Case Labels 3 - 

Explanations 

L-l/2, M-l/2 

Nominal deployment conditions; Initial orbits 
from Table 3, error magnitudes from 
Section 3, impact speed = 150 m/sec, impact 
flight path angle = -90°, AVA at apoapsis of 
initial orbit, AVB at RP or R (altitudes given 
in Table 2). 

L-3/4, M-3/4 

Sensitivities to location of AVA point; True 
anomaly of AVA reset to 170 Q from periapsis. 

M-5/6 

Sensitivities to size of initial orbit; Initial 
apoapsis altitude reduced from 24, 880 km to 
12, 140 km corresponding to a 12-hour period 
about Mercury. 

L-7/8, M-7/8 

Sensitivities to magnitude of error sources; 

O. D. and execution errors reduced by a factor 
of 3. 

L-9/10, M-9/10 

Sensitivities to nominal impact flight path 
angle; set terminal flight path angle to -80°. 

L-ll/12, M-ll/12 

Sensitivities to impact speed; increase 
nominal impact speed to 300 m/sec. 


a. Lunar cases are identified with "L", Mercury with "M"; two 

subcases are run within each case class, the first (odd number) for 
the IET Mode, and the second (even number) for the RET Mode. 
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Table 10b 

LUNAR ERROR RESULTS OF DETAILED ERROR ANALYSIS 


Cases and 
Parameters 

O.D. 
Errors 
at AVA 

AVA 

Execution 

Errors 

O.D. 
Errors 
at AVB 

AVB 

Execution 

Error 

Total 

Parameter 

Error 

L-l/2 

Rp/R(km) 

v(m/sec) 

0/12 

1/4 

— - 

— 

1/12 

-/12 

-/5 

1/- 

8/1 

8/13 

Y (deg) 

-/2 

-/14 

o/- 

15/5 

15/15 

Miss Distance, 

D m (km) 

-/14 

-/18 

15/- 

4/1 

15/23 

L-3/4 

R /R(km) 

11/5 

15/3 

— 

— 

16/6 

v(m/sec) 

-/12 

-/7 

20/- 

8/8 

22/16 

Y (deg) 

”/2 

-/13 

1/- 

15/5 

15/14 

D m (km) 

-/14 

-/IT 

14/- 

4/- 

15/22 

L-7/8 

Rp/R(km) 

o/i 

0/4 

— 


1/4 

v(m/sec) 

”/4 

-/o 

0/- 

1/3 

1/5 

Y (deg) 

-/I 

-/5 

o/- 

5/2 

5/5 

D m (km) 

-/5 

-/6 

5/- 

1/0 

5/8 

L-9/10 

Rj,/R(km) 

v(m/sec) 

0/4 

1/0 

— — 

— — 

1/4 

-/ll 

-/5 

1/- 

8/1 

8/12 

Y (deg) 

-/I 

-A3 

1/- 

15/6 

15/14 

D m (km) 

-/14 

-/18 

15/- 

4/1 

15/23 

L-ll/12 

Rp/R(km) 

0/0 

1/0 

— — 



I/O 

v(m/sec) 

-/ll 

-/5 

1/- 

4/1 

4/12 

Y(deg) 

-/I 

-A 

o/- 

8/2 

8/7 

D m (km) 

-/14 

-n 

15/- 

7/2 

16/16 
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Table 10c 

MERCURY ERROR RESULTS OF DETAILED ERROR ANALYSIS 


Cases and 
Parameters 

O.D. 
Errors 
at AVA 

AVA 

Execution 

Errors 

O.D. 
Errors 
at AVB 

AVB 

Execution 

Error 

Total 

Parameter 

Error 

M-l/2 

Rp/R(km) 

0/328 

8/14 

— 

— 

8/329 

v(m/sec) 

-/1320 

-/171 

4/- 

40/40 

41/1330 

Y(deg) 

-A3 

-/33 

o/- 

36/38 

36/52 

Miss Distance, 

D m (km) 

-/19 

-/99 

21/- 

4/4 

21/101 

M-3/4 

Rp/R(km) 

33/434 

2/646 

— 


33/778 

v(m/sec) 

-/456 

-/680 

57/- 

40/59 

70/821 

Y (deg) 

-/15 

-/30 

2/- 

36/36 

36/51 

D m (km) 

-/28 

-/92 

20/- 

4/1 

21/97 

M-5/6 

Rp/R(km) 

0/130 

8/6 

— 

— — 

8/130 

v(m/sec) 

-/144 

-/6 

4/- 

37/57 

37/155 

Y (deg) 

~/12 

-/22 

0/- 

34/36 

34/44 

D m (km) 

-/19 

-/97 

21/- 

4/4 

21/99 

M-7/8 

Rp/R(km) 

0/110 

3/2 

— 

— 

3/110 

v(m/sec) 

- / 1 1 5 

-/2 

1/- 

5/20 

5/117 

Y (deg) 

~/4 

-/H 

o/- 

12/13 

12/17 

D m (km) 

-/6 

-/33 

7/- 

1/1 

7/34 

M-9/10 

Rp/R(km) 

0/329 

8/15 

— - 

— — 

8/329 

v(m/sec) 

-/338 

-/20 

4/- 

40/13 

40/339 

Y(deg) 

-/14 

-/30 

2/- 

35/36 

35/50 

D m (km) 

-/20 

-/100 

21/- 

4/4 

21/102 

M-ll/12 

Rp/R(km) 

0/328 

8/14 

— 

— 

8/328 

v(m/sec) 

-/386 

-/17 

4/- 

20/61 

20/391 

Y(deg) 

-/7 

~/16 

o/- 

18/19 

18/26 

D m (km) 

-/18 

-/96 

21/- 

7/8 

22/98 
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R, are quite acceptable for impact speed and miss distance and marginally 
acceptable for flight path angle. Changing the nominal impact y to -80°, 
i.e. the final trajectory segment no longer being vertical, results in little 
change in final errors, except that the error in y is increased for the RET 
Mode. Increasing the nominal impact speed to 300 m/sec reduces the 
final y errors by approximately a factor of two. This is intuitively 
reasonable since the flight path angle error is primarily determined by the 
tangential velocity residual after AVB. This residual is about the same for 
both the 150 m/sec and 300 m/sec impact speed cases, but the larger 
radial component of velocity at impact with nearly the same tangential 
component reduces the flight path angle error. For the higher impact 
speed, the final y errors are close to allowable impact angles of attack, 
particularly for softer surface materials. One cautionary note needs to 
be injected at this point. The reduction of impact angle of attack errors 
by a doubling of the impact speed may be questioned, but available 
empirical data indicate that the limiting angle of attack can be expressed 
as a ratio of tangential to radial speed at impact. However, no data exists 
at impact speeds of 300 m/sec to indicate whether the ratio valid at 
150 m/sec is also valid at 300 m/sec. Indeed it may be the case that the 
maximum angle of attack may decrease with increasing impact speed in 
which case the higher speed will not alleviate the problem. 

For the lunar cases (Table 10b), similar changes in initial 
conditions produce comparable changes in the results. However, the 
initial low altitude circular orbit makes the basic situation much more 
favorable for open-loop penetrator deployment at the moon. The nominal 
cases, L-l/2, result in final errors in all parameters that are acceptable 
or nearly so. Of the input changes made, going to the higher impact 
speed, case L -11/12, resulted in a quite favorable set of final errors. 


33 



4. 4. Open-Loop Deployment Summary 

Open-loop deployment control of lunar penetrators 
produces impact errors which are very nearly within the constraints 
(Table 1). If impact soil properties are known to be fairly soft 
(S numbers > 10) and/or somewhat higher impact speeds are acceptable 
an open-loop control system would work. Conversely, open-loop deploy- 
ment at Mercury is marginal, at best. Even if an impact speed of 
300 m/sec is allowed, the IET Mode still has impact flight path angle 
uncertainties greater than 11° (the limiting impact angle of attack cited 
in Table 1, page 2 . The RET Mode is not feasible, due both to large 
impact speed errors as well as flight path angle uncertainties. 
Consequently, successful open-loop penetrator deployment is unlikely via 
either deployment mode at Mercury. The penetrator misalignment, 
between the velocity vector and penetrator longitudinal axis, is dominated 
by uncertainties in prediction of the final velocity direction. These errors 
could be slashed with onboard measurements to refine the prediction of the 
direction of the velocity vector at impact and an active attitude control 
system to realign the penetrator axis with the new predicted direction. 

The next section discusses the type of measurements and closed-loop 
control needed. 
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5 . 


ONBOARD PENETRATOR MEASUREMENTS FOR 
CLOSED-LOOP CONTROL 


The previous section on open-loop error analysis pointed to the 
need of onboard instrumentation to reduce the errors in penetrator 
deployment, particularly for Mercury missions. Two types of measure- 
ments are needed: 1) altimeter measurements to determine topographic 
relative altitudes for penetrator deployment at the AVB point, and 2) 
residual velocity measurements after the AVB maneuver to ascertain the 
direction of the impact velocity vector in order to reduce impact angle of 
attack errors. 

A radar altimeter is an obvious choice for altimetry measure- 
ments. It would be required to have a maximum range of only a few tens 
of kilometers and a minimum range of about 1 km. Range measurement 
accuracies to the order of 1% are sufficient. Small weight and power 
requirements are essential. A radar altimeter design that is quite 
acceptable in meeting these constraints is described in reference 5. A 
summary of the pertinent characteristics of this instrument is given in 
Table 11. The instrument was developed under NASA contract primarily 
for use in altimeter measurements in meteorological balloons and as such 
is not yet space qualified hardware. It could be used to give surface 
relative altitude in either the IET or RET Modes for determining the 
altitude for performing AVB. If the penetrator is released at apoapsis 
and falls to the AVB point as defined by the RET Mode, the altimeter 
would necessarily be part of the penetrator deployment system. However, 
if the spacecraft delivers the penetrator to the AVB point as defined by the 
IET Mode, the altimeter would be part of the spacecraft and could be used 
for multiple penetrator deployments. 
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Table 11 


RADAR ALTIMETER CHARACTERISTICS 21 


Mass (including antenna) 160 g 

Maximum Range >18 km 

Absolute Accuracy 1 10 m 

Power Required 0. 7 w 

RF Frequency 415 MHz 

Operating Temperature Range -55° to 55°C 


a. Reference 5 


Table 12 

PENETRATOR IMPACT ERRORS WITH RADAR ALTIMETRY 


Case and 
Parameters 

O.D. 
Error 
at AVA 

AVA 

Execution 

Error 

O.D. 
Error 
at AVB 

AVB 

Execution 

Error 

Error 

L-l/2 

v(m/sec) 

— — 

-/5 

I/- 

8/1 

8/6 

Y (deg) 

-/2 

-/14 

o/- 

15/5 

15/15 

D m (km) 

-/14 

-/18 

15 A 

4/1 

15/23 

M-l/2 

v(m/sec) 

— 

-/32 

4A 

41/41 

41/66 

Y (deg) 

-/13 

-/32 

oA 

23/38 

36/52 

D m (km) 

-/19 

-/99 

21 A 

4/4 

21/101 
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Table 12 gives the effect of including altitude measurements on 
the terminal errors. The altitude measurements were simulated by 
reducing a of the AVB altitude to 0. 03 km in the O. D. error matrix and 
in the error covariance resulting from propagation of AVA errors to the 
AVB point. Correlations between position and velocity must also be 
considered. For example, measurement of altitude means that the 
position of the penetrator is determined to an accuracy limited primarily 
by uncertainties in planetary topography. Knowing the altitude of the 
penetrators permits determination of its speed to an accuracy consistent 
with initial error sources. For Mercury, the speed prior to AVB is 
insensitive to the assumed O. D. error sources and to errors in Mercury's 
topography on the order of 1 km. However, execution of AVA introduces 
a tangential velocity error which propagates to nearly 103 m/sec at the 
AVB point. This error, together with AVB execution errors, is now the 
dominant source of impact errors. It can be seen in Table 12 that the 
addition of altitude measurements results in a substantial reduction in 
impact speed errors but has little effect on other errors. This is 
because impact speed errors, particularly in the RET Mode, result from 
errors in the altitude of the AVB point which are significantly reduced by 
making altitude measurements. Errors in flight path angle, Y, and miss 
distance, D^, are dominated by errors in the horizontal velocity. These 
errors are essentially unaffected by altitude measurements. 

Determination of the flight path angle, y, at impact requires 
knowledge of the relative horizontal ground speed, v^., of the penetrator 
after AVB has been performed. If v^. is measured immediately after 
AVB, which occurs at an altitude H, then y at impact is determined from: 
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tan y = 


(3) 


* ^ > 

/i H x 

(i + H )v t 

S 

where v is the radial velocity component and R is the surface radius of 
rs s 

the planet. Using the nominal value of v in equation (3), then allows y 

r s 

to be predicted. The horizontal speed may be measured in one of two 
ways: 1) onboard measurement of the maneuver execution velocity errors, 
or 2) direct measurements of the horizontal speed relative to the surface. 

5. 1. Closed-Loop Attitude Control with Accelerometers 

The first method is possible since the largest contribution 
to the residual horizontal velocity arises from execution errors in per- 
forming AVB, hence measurement of the velocity residuals gives a good 
estimate of v^.. These errors could be measured through the use of 
onboard accelerometers. To ascertain measurement accuracies of 
onboard accelerometers for penetrator deployment the representative 
error characteristics given in Table 13 were used. 

Calculation of velocity measurement accuracies also requires 
knowledge of the accelerations involved in performing the maneuvers. 
Several solid propulsion systems having mass of the order required for 
AVB at Mercury have been developed for space applications. Thrust levels 
on representative systems vary from about 17, 000 newtons to 170, 000 
newtons, with burn times from around 3-20 seconds and would give peak 
accelerations of 33-330 g's. Table 14 gives the resulting errors in 
velocity measurements of AVB of the IET Mode at Mercury for the high 
thrust level (170, 000 newtons) based on the error sources of Table 13. 

The smaller thrust level would result in increasing the bias error to 
3x10” m/sec, which is still negligibly small. Although various types of 
accelerometers have different error characteristics, the above analysis 
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Table 13 

TYPICAL ACCELEROMETER ERROR CHARACTERISTICS 

-4 

Bias 1x10 g’s 

-5 

Scale Factor 5x10 g’s/g 

Misalignment 80 arc sec 

Table 14 

ACCELEROMETER MEASUREMENT ERRORS OF AVB a 



Bias 

Scale 

Factor 

Misalignment 

Ax(m/sec) 

0 

0 

1.6 

Ay(m/sec) 

3xl0“ 3 

0.2 

0 

Az(m/sec) 

0 

0 

1.6 


a. Mercury IET Mode; nominal AVB = 4065 m/sec 


Table 15 

IMPACT ANGLES OF ATTACK USING ACCELEROMETERS a 


Case and 
Parameter 

O.D. 
Error 
at AVA 

AVA 

Execution 

Error 

O.D. 
Error 
at AVB 

AVB 

Execution 

Error 

Total 

Error 

Ml 

a(deg) b 

. 


0.09 

1.12 

1.12 

LI 

a(deg) 

- 

- 

0 

1.10 

1.10 


a. IET Mode only 

b. a = impact angle of attack 
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indicates that the velocity measurement errors will be at most 
2-3 m/sec per axis. Measurements of this accuracy are quite sufficient 
to reduce the effects of maneuver execution errors on the impact angle of 
attack to acceptable levels. 

In Section 4 it was noted that there are two components to the 
residual tangential speed: 1) an in-plane component in the direction of the 
AVB vector arising from impulse errors in the burn magnitude, and 
2) an out-of-plane component due to thrust misalignment. These 
residuals are of nearly equal magnitude and both would have to be 
measured. This requires two accelerometers and probably a three 
accelerometer system would be simplest to mechanize with gyros. A 
strapdown system with the accelerometers and gyros attached directly 
to the penetrator frame would be well suited as this system offers weight 
saving at the expense of some measurement accuracy. 

Table 15 summarizes the impact angle of attack errors for the 
IET Mode based on accelerometer measurements of all three velocity 
components and the measurement accuracies described above. Following 
the AVB burn, the penetrator reorients to the path angle and azimuth 
determined by measuring velocity components. Assuming 3a final 
attitude control pointing errors of 1°, the 3a angle of attack error at 
either the moon or Mercury would be less than 2°. This error is well 
within the 3° constraint for impacts in very hard soils, e.g. basaltic 
lava. Note, however, that the impact flight path angle errors presented 
in Table 12 are unchanged by either the radar altimeter or accelerometer 
measurements. While penetration can probably be achieved, the resulting 
penetrator angle at rest may be unacceptable for some experiments, 
particularly at Mercury where the 3a y errors are considerably larger 
than the desired 15° limit cited in Table 1, page 2 . 
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5.2. 


Closed-Loop Attitude Control with Driftmeters 


The second method of measuring the residual horizontal 

velocity after AVB is a passive optical technique employed by driftmeters 

(v^/H meters). Specifically, optical driftmeters measure the ratio v^/H, 

where v is the horizontal velocity at the altitude H. These devices have 

also been considered for determination of attitude rates of orbiting 
7 

spacecraft . The basic technique consists of an optical system to form 
the image of the surface onto a reticle consisting of alternate opaque and 
transparent bands and then onto a photodetector. Motion of the v^/H 
meter parallel to the surface moves the image across the banded surface, 
which effectively chops the image to produce an output signal whose 
average frequency is proportional to v^./H. Measurements of v^/H can 
then be combined with radar altimeter measurements of H to determine 
v t and consequently allow prediction of y. 

The measurement accuracy required in order to meet 
the most stringent angle of attack constraint of 3° may be estimated by 
differentiating equation (3) above, yielding 


6v 


sec y 5>y = 


rs 


(1+ 1 )v t 
s 


v 6v, 
rs t 

H v 2 

(1+ R )v t 
s 


V 

rs 


n H .2 
5 S v ‘ 



( 4 ) 


For Mercury, H/R =. 0012 and 5— » 5x10 , assuming for 6 H, a 

s Kg 

conservative estimate of the 3 ct altimeter error of 0. 1 km. Neglecting 
these terms, equation (4) simplifies to 
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( 5 ) 


2 6v 

sec y5y w rs 


't 


v 

rs 



For these terms to be approximately equal in magnitude requires that 

Sv t _ ^ V rs . Using a value of 150 m/sec for v and an uncertainty of 

v, v 
t rs 

18 m/sec for 5>v results in ^ v t ~0. 12. These values come from the 
rs 

V t 

nominal Mercury IET Mode where the 18 m/sec error in radial impact 
speed results from a 3a residual radial velocity of 71 m/sec after AVB 
(see Table 7, page 20). The uncertainty in the 6v^. measurement arises 
from errors in measurement of v^/H, and 6(v^./H), together with altitude 
measurement uncertainties 6H. Since 


6v, 6(v./H) 6H 

^- = 57H— + tt’ 


and our estimate of 6H/H is 0. 03 at Mercury, then for a value of 0. 12 for 
6v^/vj. equation (6) provides 


6(v./H) 


0.09. 


Hence, a driftmeter measurement accuracy of about 10% is sufficient to 
reduce the uncertainties in the predicted value of yto the same level as 
those resulting from the unmeasured radial speed residual after AVB. 

The above discussion is valid only if the residual tangential speed 
after AVB is aligned with the penetrator mounted driftmeter. In general, 
errors introduced by the AVB burn will also contain an out-of-plane 
component in the tangential velocity, v^., perpendicular to this direction. 
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Therefore, two orthogonal driftmeters are required to sense the true 
magnitude of v^.. The impact attitude (angle of attack) errors resulting 
from closed-loop attitude control using a pair of driftmeters are 
summarized in Table 16 for the IET Mode at the moon and Mercury for 
three levels of driftmeter accuracy: 5, 10, and 15%. Note that these 
results are based on a nominal impact velocity of 150 m/sec, a 3a error 
in the radial impact speed, and an inherent 1° error in the attitude control 
system pointing accuracy. Results are given for three levels of residual 
tangential velocity. From Table 7, page 20, 3a residual tangential 
velocities at the moon and Mercury are expected to be 39 m/sec and 
94 m/sec, respectively. Comparing these values with the data in Table 16, 
driftmeter supported closed-loop attitude control can be expected to 
reduce impact angle of attack errors to 4° at the moon and 9° at Mercury. 
Recalling that closed-loop attitude control with an accelerometer triad 
reduces angle of attack errors to about 2° at either target, it can be seen 
that accelerometers are preferred to driftmeters, particularly if harder 
impact soils are expected where the maximum angle of attack error 
should be less than 3°. This conclusion is strengthened by the fact that 
accelerometers are completely independent of target surface features 
(shadowing) and provide results which are less sensitive to the radar 
altitude measurement. Note that the angle of attack errors using 
accelerometers is nearly the same as the systematic error introduced in 
Section 2 by ignoring the rotational velocity of the primary. However, 
targeting the final trajectory segment to allow for the finite rotational 
speed would eliminate the systematic error but there would still remain 
the random angle of attack errors. 
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Table 16 


IMPACT ANGLES OF ATTACK USING DRIFTMETERS a 


6(v./H) 

Target ~v^7h~ 

v t 

2v t 

3 v t 

• Moon (IET Mode, v^. = 13 m/sec) 




0.05 

1.2° 

1.7° 

2,3° 

0.10 

1.6° 

2.6° 

3.7° 

0.15 

2.1° 

3.7° 

5.3° 

• Mercury (IET Mode, = 30 m/sec) 



0.05 

3.0° 

5.3° 

6.8° 

0.10 

3.8° 

6.9° 

8.9° 

0.15 

4.8° 

8.8° 

11.5° 


a. Results include a 1° attitude control pointing error 
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6 . 


SUMMARY AND CONCLUSIONS 


The error analysis of penetrator deployment from orbiting space- 
craft at targets without atmospheres yields a number of significant results. 
First, the Intermediate Ellipse Transfer (IET) Mode is preferable both 
from the viewpoint of minimum AV requirements and deployment errors. 
Second, open-loop error analysis indicates that when the penetrator is 
deployed at apoapsis and coasts on its own to the AVB maneuver point, 
large errors occur in the impact speed and flight path angle. At Mercury 
the errors are unacceptably large, whereas for lunar deployment the 
flight path angle errors are marginally acceptable for impacts into soft 
material, but are too large for impacts into rocky surfaces. The 
terminal errors in impact speed are significantly reduced if the space- 
craft can deliver the penetrator to the second maneuver point, thereby 
eliminating the effects of initial O. D. errors and maneuver execution 
errors at AVA. For Mercury orbiters it may be possible to utilize solar 
perturbations on the orbit to smoothly reduce periapse at small rates in 
order to attain the low 3 km periapse altitude required for penetrator 
deployment. While this deployment mode is preferred from the viewpoint 
of reducing penetrator errors, the effect of orbit orientation constraints 
on the overall mission design must be examined. The largest open-loop 
error is in the angle of attack at impact. This error arises primarily 
from errors in the predicted flight path angle at impact due to the 
horizontal velocity errors introduced by the AVB maneuver. These errors 
are due to out-of-plane pointing errors in performing AVB and to errors 
in the magnitude of the AVB impulse. A second problem area for Mercury 
deployment is the magnitude of the flight path angle at impact. This angle 
should be no greater than 15° from the vertical, whereas the 3a errors are 
36°. This can only be controlled, in the absence of propulsive maneuvers, 
by using a higher impact speed or by reducing the impulse variations of 
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the rocket engine and out-of-plane pointing errors. Finally, the errors 
in impact speed are near the limits of acceptability. These errors can be 
reduced only by reducing the error sources in AVB just noted, again 
assuming no corrective propulsive maneuvers. 

Orbiter altitude measurements are required for precise deter- 
mination of the AVB maneuver altitude. For both lunar and Mercury 
penetrator deployment an onboard determination of the impact flight path 
angle is required in order to reduce penetrator angle of impact attack 
errors. Onboard accelerometers measuring the velocity residuals after 
AVB provide an accurate measurement of the primary error source for the 
impact angle of attack. An alternative onboard measurement technique 
using optical driftmeters, in conjunction with altimeter measurements, 
cannot reduce the angle of attack errors to acceptable levels for all cases. 

Although higher impact speeds are desirable from the viewpoint of 
reducing errors in the impact flight path angle and allowing higher rest 
altitudes, they quickly lead to difficulties with the penetrator design. As 
discussed in Section 2, doubling the impact speed from the nominal value 
of 150 m/sec to 300 m/sec would exceed the deceleration limit of 2000 g’s 
in rocky soil and the penetration depth of 15-20 m in soft soils. Changing 
the penetrator diameter would solve the problem at one end but make it 
worse at the other end. Also, the penetrator forebody diameter of the 
Mars design is only 9 cm, so making it any smaller to ease g-loads in 
hard soils would also create new payload packaging problems. 

One final area of concern was identified but not analyzed. It 
involves the response time of the attitude control system. Immediately 
after the AVB maneuver the onboard control system of the penetrator 
must reduce the accelerometer measurements to the desired penetrator 
orientation, nominally about 90° from the horizontal AVB maneuver 
attitude. Then the attitude control system must pitch the penetrator to 


46 



this attitude, damp out over /under- shoot errors and do so with enough 
time to jettison all systems external to the penetrator itself (31 kg) before 
impact. This jettison must occur sufficiently early and with enough speed 
to guarantee that impacts of the excess parts do not interfere with the 
penetrator afterbody or its antenna, either by hitting it directly or through 
ejecta created by their impacts. The nominal fall times for a 150 m/sec 
impact speed, are 41 sec at Mercury and 93 seconds at the moon. Using 
the 3a estimates of the error in the radial component of velocity from 
Tables 6 and 7, (page 20), and assuming that the error is directed toward 
the surface, the fall times drop to 26 seconds at Mercury and 76 seconds 
at the moon. Whether or not this is sufficient time to accomplish all the 
operations just described successfully is a question of feasibility which 
has not been analyzed. To do so will require a more detailed definition 
of the penetrator design and its mass properties for Mercury missions. 

In conclusion, the deployment of lunar penetrators appears to 
pose no unreasonable performance or control requirements. Conversely, 
the low deployment altitude, large retro mass, large retro execution 
errors, and short descent time, all raise fundamental feasibility questions 
for a Mercury penetrator mission. For either concept more detailed 
design and error analysis will be required before serious mission 
planning can reasonably be considered and is recommended. 
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APPENDIX A: FORMULATION OF PARTIAL DERIVATIVES 


To determine the effects of the various error sources on the final 

variables of interest (periapse radius, radius magnitude, speed, flight 

path angle and miss distance), the partial derivatives of the final 

variables with respect to the initial position and velocity components are 

needed. Fundamental to the calculation of several of these derivatives is 

the state transition matrix, which is a 6x6 matrix of the partial derivatives 

of the position vector, r, and velocity vector, v, with respect to an initial 

position vector, r , and initial velocity vector, v q , along a Keplerian 

trajectory. A formulation of the state transition matrix in universal 

variables, based upon reference 6, was utilized to compute the partial 

derivatives, “ = ~ , etc.^ , etc. for all of the 

’ Zr lax ’ ay ’ dz ’ Zx * } ’ 

required 36 elements. These derivatives were available for use in the 

calculation of the partials derived below. 

Periapse Partial Derivatives 


The derivatives of periapse radius, RP, with respect to the initial 
state vector, • - 3RP aHP 3RP 3RP 3RP -- 3RP 


i. e. 


dx. 


ay. 


dz. 


required. 


o o " o 
Using RP = a (1-e), then 


0x 


o 


dy r 


and 


Zz 


were 


o 


^RP 

^ X o 


aa 

ax 

o 


(1-e) 



A. 1 


From the energy relation 


2 

v p u . 

2 r "r- -2?> onehas 


A. 2 


o 2 

da _ 2a xo . 
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The term may conveniently be evaluated in terms of the partial 

ZXq 

derivatives of two other functions, C and S , which are tabulated in 

’ e e’ 

reference 8, p. 331. Now from the definitions C g = e cos E, S g = e sin E, 
where E is the eccentric anomaly, we have 


3x o 


ae 

ax o 


cos E -e sin E 


aE 

3x o 


A. 3 


as 


e de . „ _ aE 

— - — sin E + e cos E 


ax„ ax^ 

o o 


ax. 


A. 4 


Multiplying A. 3 by cos E and A. 4 by sin E and adding gives 


ae 

3x o 


3C e 

cos E + 

ax Q - o 


as 

----- sin E 
ax o 

o 


From reference 8, 



x v 
o o 


9 


as x 

e_ _ o 

ax i — 

o Vua 


aS x 
e o 



Combining A. 1, A. 2, A. 5, A. 6 and A. 7 yields 


A. 5 


A. 6 


A. 7 
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dRP 


2a x (1-e) 
o 


dx 


- a 


o 


o 


x v cos E + 
o o o 


\ir 


o 


/ * 


sin E 


o 


aS x 
e o 


& 


A. 8 


The partials , may be found from A. 8 by replacing 

y o z o 

(x q , X Q ) by (y o , y Q ) and (z q , z q ) successively. 

The velocity partials may be found in an anologous manner. The 
energy relation may again be employed to determine 


o 2 - 

2a x 

da _ o 


dx Q P 


Then 


dRP 

dx„ 


2a^x (1-e) 
o 

P 


-a 


2r x cos E 
o o o 


P 


/x. 


sin E 


VET 


aS x 
e o 

P 


A. 9 


Again the remaining partials, — — — and — are found by 


replacing (x q , x q ) with (y Q , y Q ) and z q , z q ). 
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Radius and Velocity Partial Derivatives 


2 2 2 2 

From the relation R = x + y + z 


R 


x o R 


hx by &x 

X -rr— + y + Z 


dx 


As the elements 


o 

~dx 

dx 


hx 


o 


dx 


A. 10 


, etc. are available from the state transition 


matrix, these partial derivatives of R are straightforward. Similarly, 
using 

2 .2 .2 .2 

v = x + y + z , 


then 


dv 


dx 


o 


1 

V 


. dx 

X 




+ y - 


^X, 


. b z 

+ z 


ax. 


A. 11 


and the derivatives , etc. are obtained from the state transition 

■v " 


X 

o 

matrix. The other partial derivatives are obtained by replacing x q by 
y Q , z q , x q , y Q , and z q successively. 

Flight Path Angle Partial Derivatives 

From the definition of the flight path angle y , one has 


sin y = 


r . v 
rv 


Forming the partial derivative with respect to x q , 


cos y 


'by _ 1 ^(r . v) r . v ^(rv) 


ax„ rv 3x, 


(rv) 


2 ax. 


A. 12 
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Now 


j-Kg - y) = ^>(xx + yy + zz) = . ^>x ^x . V 

ax o *x o ax o ; S x o y ix o 


'hy • ^ z ^z 

+ y — - — + Z -r-r— + Z 




o 


ax f 


2>x, 


A. 13 


Equation A. 13 may be evaluated using the state transition matrix. The 
term in Equation A. 12 can be expressed as 


Mrv) = v + r 


* x . 


dX, 


* X o 


A. 14 


Note that ^ r and have been previously evaluated in Equations 


d X , 


b x 


o v o 

A. 10 and A. 11. Hence, combining A. 12, A. 13 and A. 14 yields 


»Y _ 


2>x, 


rv cos y 


( 


x> 


x 


* x . 


+ X 


^>x 

ax. 


• 'by 
+ y 


d x 


+ y 


*y 

^X 


A. 15 


. 'bz 
+ z 


dx 

o 


+ z 


^>z 

*x. 


\ (xx+yy+zz) f br b 

j ” / n2 Jx ^ 

i / (rv) cos y ' ^ y 


y-\ 

x I 

o / 


The other partials may again be found by successively replacing 

(x , x ) in A. 15 by (y , y ) and (z , z ). 

\ o’ o J J o J o v o o 


Miss Distance Partials 


Let the polar co-ordinates of the impact point be a, 6, where 
sin b = Z /Rpi ane j- and tan a = y/x with x, y, z being planetocentric 


cartesian co-ordinates. 
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Then 


and 


cos 6 


*6 


I 


^ x o ^planet ^ X o 


„ 2 
sec a 


by/l* 


o 


h x . 


2 

x o 


or, 


"6a 2 

<- = cos 


a(i|Z- - 

\ x ^ x o x 2 ix o I 


Now for a small change Ax q , the change in the impact point or 

miss distance (D ) is: 
v m 


A D — R , . , A 8 — R •. , 

m planet planet 


[ 


(A6) 2 + cos 2 6 (Aa) 


■] 


1/2 


= R 


planet 


[■* 


v2 2 , ^ a v 

) + cos h ( ) 


2 ] 1/2 

^ X o ' J 


Hence, 






= R *9 = R 

planet ^x q planet 


, x2 2 / ^ 

{ jr- ) +cos 6 ( s 


o J 


1/2 


A. 16 


The other partials are obtained as before. 


Equations A. 15 and A. 16 become indeterminate (zero over zero) 
along rectilinear trajectories. In such cases the partial derivatives are 
computed numerically. 
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APPENDIX B: ERROR ANALYSIS COMPUTER PROGRAM 


A program was developed on the CDC 6400 computer to compute 
the maneuver AV's, O. D. and execution errors, and errors in the 
dependent variables. The required program inputs are given in Table B-l. 
The program first computes the nominal AVA and AVB maneuvers, the 
intermediate and terminal trajectories and impact conditions based on 
conic trajectories. From the position and velocity at the AVA and AVB 
points and AVA and AVB magnitudes, O. D. and execution error covariance 
matrices are calculated. The errors are assumed to be uncorrelated, 
hence these are. diagonal matrices. Next, the partial derivatives are 
evaluated along the trajectory after the AVA maneuver has been performed. 
The formulation of partial derivatives is given in Appendix A. Based on 
these partials and the O. D. and execution errors at the AVA point the 
error in the AVB maneuver radius is determined. The O.D. covariance 
and AVA execution covariance matrices are also propagated to the AVB 
point utilizing the state transition matrix evaluated along the intermediate 
trajectory. The state transition matrix is based upon a formulation in 
universal variables (and consequently is valid along any conic trajectory) 
as described in Reference 6. Finally, the program calculates the partial 
derivatives along the final trajectory segment and determines errors in 
the impact parameters. 

Program output consists of conic trajectory parameters and AV's, 
error sources for O. D. and execution errors, partial derivatives along 
both trajectory segments and resulting errors in periapse radius, impact 
speed, impact flight path angle, impact angle of attack, and miss distance. 
Table B-2 lists the specific program outputs. 
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Table B-l 


ERROR ANALYSIS PROGRAM INPUTS 


Variable Name 

Units 

Description 

GMU 

i 3 / 2 

km /sec 

Gravitational parameter of the 
target body 

RAO 

km 

Apoapse radius of initial orbit 

RPO 

km 

Periapse radius of initial orbit 

RPLAN 

km 

Radius of target body 

MODE 

- 

= 1 for I. E. Mode 
= 2 for R. E. Mode 

VINOM 

km/sec 

Nominal impact speed 

GAMINOM 

deg 

Nominal impact flight path angle 

RPINT 

km 

Periapse radius of intermediate 
ellipse for I. E. mode 

RINT 

km 

Radius of AVB maneuver point 

TAO 

deg 

True anomaly of AVA maneuver 
point in initial orbit 

EXERR(2) 

- 

Fractional error in magnitude of 
AVA and AVB 

OOPPE(2) 

deg 

Out of plane pointing error for 
AVA and AVB maneuvers. 

IPEE(2) 

deg 

Inplane pointing error for AVA and 
AVB maneuvers 


58 



Table B-l (Concluded) 

ERROR ANALYSIS PROGRAM INPUTS 


Variable Name 

Units 

Description 

ODCVCHG 

- 

Numerical multiplier used for 
changing magnitude of O. D. errors 

SIGATT 

deg 

3ct error in attitude control system 
pointing errors 

IENDFL 


Set to 1 on last case to terminate 
run 
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Table B-2 


ERROR ANALYSIS PROGRAM OUTPUTS 


A. NominalConditions 

Initial position and velocity vectors 
AVA and AVB vectors 
Position and velocity vector at AVB 
Velocity vector after AVA and AVB 

Semimajor axis, eccentricities and flight times along both 
trajectory segments 

B. Error Sources 


Orbit determination error matrix at AVA point, AVB point and 
AVA errors propagated to AVB 

Execution error matrix for AVA and AVB maneuvers and AVA 
errors propated to AVB 

C. Dependent Variable Sensitivities 

Partial derivatives of R, RP with respect to initial position and 
velocity vector components computed along the trajectory from 
AVA to AVB 

Partial derivatives of impact speed, flight path angle and miss 
distance with respect to position and velocity vector components 
along the trajectory from AVA to impact 

D. Error Independent Variables 


Total error and errors due to individual sources are computed 
for: 

o AVB maneuver altitude 
o Impact speed 
o Impact flight path angle 
o Impact angle of attack 
o Miss distance 



